Correlations of the local density of states in quasi-one-dimensional wires 
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We report a calculation of the correlation function of the local density of states in a disordered 
quasi-one-dimensional wire in the unitary symmetry class at a small energy difference. Using an 
expression from the supersymmetric sigma-model, we obtain the full dependence of the two-point 
correlation function on the distance between the points. In the limit of zero energy difference, our 
calculation reproduces the statistics of a single localized wave function. At logarithmically large 
distances of the order of the Mott scale, we obtain a reentrant behavior similar to that in strictly 
one-dimensional chains. 

PACS numbers: 73.20.Fz, 73.21.Hb, 73.22.Dj 



I. INTRODUCTION 



Anderson localization in quasi-one-dimensional (Q1D) 
disordered systems is known as a broad universality class 
for various problems of contemporary condensed mat- 
ter physics and of quantum chaos 
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luasi-one- 



dimensional wires are characterized by a large number 
N ^S> 1 of transverse channels. Starting from the mean 
free path scale I, electrons propagate diffusively and 
eventually get localized^ with the localization length 8 
£ ~ Nl 3> I. Quantum diffusion in the broad range 
of length scales is the key feature of Q1D systems dis- 
tinguishing them from strictly one-dimensional (ID) dis- 
ordered chains, where localization occurs immediately at 
the scale of the mean free path^il^ It is this intermediate 
diffusive regime that makes studying localization in Q1D 
systems a very challenging problem. 

In order to describe the nonperturbative regime of 
strong localization, a rather sophisticated mathematical 
technique is required, both in the ID and Q1D cases. For 
the strictly ID geometry, an appropriate technique has 
been elaborated by Berezinsky , 11 ' 12 and recently it was 
translated into a field-theoretic language in Ref.hJ. In a 
Q1D geometry, where localization takes place in the pro- 
cess of quantum diffusion, the diffusive supersymmetric 
cr-model introduced by Efetov^^ is a powerful tool to 
describe the system. 

Despite the physical difference between the two mod- 
els, localization in the ID and Q1D geometries often 
look similar. This analogy is most pronounced in the 
behavior of a single localized wave function: The statis- 
tics of the smooth wave function envelopes are precisely 
the same in the I D 15 : 16 ' 17 and Q1D 18 geometries, while 
the local statistics of the rapidly oscillating components 
involves the plane-wave (ID) or random-matrix (Q1D) 
moments* 1 ^ 

On the other hand, comparison between the quanti- 
ties involving different wave functions in the two models 
is problematic since almost nothing is known about the 
eigenfunction correlations in the Q1D geometry in the 



most physically interesting limit of sufficiently small en- 
ergy difference u. 

At the same time, quite a lot is known about eigen- 
function correlations in the strictly ID geometry un- 
der similar conditions. One of the most important re- 
sults is the expression for the dissipative low-frequency 
conductivity which follows the Mott-Berezinsky law i 11 ! 20 
Recr(w) oc lu 2 ln 2 (l/cjr), where r is the elastic time. An- 
other result concerns the behavior of the correlation func- 
tion 



i?(w;ri,r 2 ) =v 2 {p e (Y 1 ) Pe+u) {v 2 )) 



(1) 



of the local density of states (LDOS) Pe( r ) = 
J2 n IV'n( r )| 2( 5(e — £n), where v is the average density 
of states (for spinless particles). The dependence of 
i?(w;ri,r 2 ) on |ri — r 2 | has been addressed for ID ge- 
ometry by Gor'kov, Dorokhov and Prigara 2 --. They have 
identified that in the limit of small energy separation, 
Mr C 1, the spatial behavior of the LDOS correla- 
tor is determined by two scales: the localization length, 
£id ~ I, and the Mott scale, L 1 ^ ~ £id \n(l /u>t) S> £id- 
The eigenstates are uncorrelated at r ^ and exhibit 
nearly perfect level repulsion at £id <C r <C Lff, which 
gets reduced at r <C £id- Such a behavior can be qual- 
itatively understood 22 using the Mott's picture^ 3 , of two 
nearly degenerate localized states separated by the dis- 
tance Ljtf which hybridize to form a pair of eigenstates 
effectively contributing to R{w; r\, v%). 

In the Q1D geometry, neither the derivation of the 
low- frequency conductivity nor the full dependence of the 
LDOS correlation function have been reported. Though 
the one-dimensional diffusive supersymmetric cr-model 
has been reformulated by Efetov and Larking as a 
quantum-mechanical problem, the resulting set of differ- 
ential equations is still too complicated and resisting a 
naive perturbative expansion in the limit to — ► 0. 

Recently, a significant progress in solving this effec- 
tive quantum mechanics has been achieved by two of the 
authors^! who found the exact zero mode of the transfer- 
matrix Hamiltonian for the unitary symmetry class. As 



2 



a result, a non-perturbative expression for the short-scale 
(at distances |r x — r 2 | <C £) behavior of the LDOS cor- 
relator JT]) has been derived for arbitrary frequencies u>. 
For the first time, the difference in the correlations of the 
localized states in the ID and Q1D geometries has been 
demonstrated. 

The purpose of the present work is to develop a regular 
approach to the small-w expansion of the one-dimensional 
diffusive supersymmetric er-model. We will assume the 
unitary symmetry class and employ the exact knowledge 
of the zero mode obtained in Ref. [27|. Developing an ad- 
vanced perturbation theory which treats both powers and 
logarithms of u, we will derive the expression [Eq. ([66)) ] 
for the LDOS correlation function at arbitrary distances 
|ri — r2 1 between the observation points and analyze its 
behavior in various asymptotic regions. 

The LDOS correlation function (TTJ) can generally be 
written as2£ 




FIG. 1: (color online). Dependence of the LDOS correlation 
function R(u>, x)ont = x/£ for various values of w/Aj shown 
on the graph. 



fl(w;ri,r 2 ) = 1 + A(u,t) + k(r u r 2 )B(uj). (2) 

Here the factor fc(ri,r 2 ) = (ImG fl (ri,r 2 )) 2 /(7r^) 2 ac- 
counts for short-scale Friedel oscillations ! 19 ' 31 It is equal 
to 1 at coincident points and exponentially decays as 
| r i — r*2 1 exceeds the mean free path. The function A(u>, t) 
describes long-range correlations, with 



(3) 



measuring the distance x — x\ — x 2 between the obser- 
vation points in units of the localization length, 



£ = 27WiL>, 



(4) 



where D is the diffusion coefficient and v\ = vS is the 
one-dimensional density of states (S is the wire cross sec- 
tion) . 

The functions A(u>, 0) and B(oj) which determine local 
correlations at x -C £ have been calculated in Ref. 



A{u, 0) = - Re [n 2 (I 2 - I I 2 ) (K 2 - K K 2 ) - I 2 Kf] , 



B(u) = -Re(l 1 K 1 +2I 2 K ), 



(5) 
(6) 



where k = yf— 4iw/A^ and the argument k of the modi- 
fied Bessel functions is suppressed. 

At large separations exceeding the mean free path, x 3> 
I, the function fc(r l7 r 2 ) exponentially decays and Eq. @ 
simplifies to 



R(w,x) = 1 + A(u,t). 



(7) 



The behavior of R(ui,x) as a function of t — x/£ cal- 
culated with the help of the general expression ([55)1 is 
shown in Fig. [T] Different curves on the graph corre- 
spond to different values of the ratio lo/A^, where the 
level spacing within the localization length, 



a, = Die 



(8) 



is the natural frequency scale in the localization problem. 

One can clearly see that in the deeply localized regime, 
uj -c , the behavior of R(ui, x) is similar to the behav- 
ior in the strictly ID geometry; 21 ' 22 On increasing x, the 
function R(u>,x) first decays at the localization length 
and then reaches the uncorrelated value R = 1 at the 
Mott scale 



L M = 2£ln(A £ /u;). 



(9) 



The function R{lo, x) has a minimum of the order of 
(w/A^) 1 / 4 at x~ L M /2. 

The paper is organized as follows. In Sec. [Til starting 
with the supersymmetric er-model formalism, we give the 
expression of the LDOS correlation function in terms of a 
matrix element of the evolution operator with the effec- 
tive Hamiltonian. The bosonic and fermionic eigensys- 
tems of this Hamiltonian are perturbatively constructed 
in Sees. [HTJ and IIV1 respectively. In Sec. [V] we give the 
general expression for the LDOS correlator and evaluate 
it for arbitrary x and small to. The opposite limit of 
small x and arbitrary lo is considered in Sec. IVIl Section 
IVIII is devoted to the discussion of results and compar- 
ison between the ID and Q1D localization. Finally, a 
large portion of technical details is relegated to several 
Appendices. 



II. SUPERSYMMETRIC SIGMA-MODEL AND 
GENERAL EXPRESSION FOR THE LDOS 
CORRELATION FUNCTION 

An efficient analytical tool for studying disordered 
metallic systems — nonlinear supersymmetric er-model 
- was proposed in the works of Efetovii 4 -. This ap- 
proach is based on the effective action 



S[Q] =-^-str J dx 



D(^) r 2/,,AO 



(10) 
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for slow diffusion modes described by the 4x4 super- 
matrix Q acting in the direct product of the Fermi-Bose 
(FB) and Retarded-Advanced (RA) spaces (we assume 
the unitary symmetry class). The Q matrix is sub- 
ject to the nonlinear constraint Q 2 — 1. In Eq. (p~0|) . 
A = diag(l, — I) the matrix in the RA space. 

The one-dimensional cr-model (fTUf can be mapped 2 ^ 
onto a quantum mechanics of the Q matrix, with the 
coordinate x along the wire playing the role of time. In 
that language, the functional integral with the action (fTU]) 
describes evolution of the wave function ^f(Q) under the 
action of an appropriate Hamiltonian H,: 



|Ukq) = -2w*(q), 



(11) 



where t is the dimensionless "time" ([3]) measured in units 
of the localization length. In the unitary symmetry class, 
the Q matrix is generally parameterized by four com- 
muting and four Grassmann variables. However in many 
cases it suffices to consider only the "radially-symmetric" 
("singlet") wave functions ^(Q) = ^(X B ,X F ) which 
depend only on Efetov's coordinates Xb G [L°o) and 
X F G [—1,1], parameterizing the non-compact BB and 
compact FF sectors.*!*^ On the space of such functions, 
the Hamiltonian Ti acquires the formi^i 



H 



(Ai 



X F ) 2 



d 



1 



d 



d 



dX F {X B - X F ) 2 dX F 
- 1 d 



+ n(x B -x F ), (12) 



dX B (X B - X F ) 2 dX B _ 
where we have introduced the dimensionless frequency 



4A 



(13) 



At u> = 0, the Hamiltonian (| 1 2[) reduces to the Laplace 
operator on the space of Q-matrices. This operator is 
know n 25 ' 26 to have a discrete spectrum in the FF sector 
and a continuous spectrum in the BB sector. The main 
technical difficulty in treating the finite-cj case is that 
any nonzero to renders the continuous BB spectrum be 
discrete (assuming 51 positive). This is related to the 
noncompactness of the BB sector eventually responsible 
for localization. 

The zero mode of the Hamiltonian (fT2|) solving H^o = 
has been recently obtained in a closed form valid for all 
frequencies^ 



*o(Af, X b ) = h{q)pKi[p) + qh(q)K (p), 



(14) 



where 



P = y/8si(\ B + i), 



vmxi+T). (15) 



The function A(u>,t) which determines the large- 
distance behavior of the LDOS correlation function ([7]) 
can be conveniently written in terms of the expectation 



value of a certain evolution operatori 2 ^ calculated with 
respect to the zero mode (fT4"|): 



A(u,t) = iRe<* |e- 2 ^|v]/ ), 
where (-|-) is the fiat bilinear form on singlet states: 



(16) 



(<f>\x) 



dX f 



dXs <)>(X B , X f )x(Xb, X f ). (17) 



The Hamiltonian H = (Xb — X F ) 1 H(Xb — X F ) splits 
into the sum of the bosonic and fermionic parts: 



H = H F 



H F , 



where 



H B = 
H F = 



-^d XB (X%-l)d Xl 



nx f 



-^d XF (i-x 2 F )dx F -nx F . 



(18) 

(19) 
(20) 



Analytic calculation of the expectation value in (QI 
as a function of arbitrary lo and t is a complicated task. 
The main technical problem is that the eigenfunctions of 
the Hamiltonians Hb and H F (known in the theory as 
a particular class of the Coulomb spheroidal functions 29 ) 
cannot be found explicitly for an arbitrary value of the 
parameter O. 

In the limit Q 3> 1, only small deviations of and 
Xb from 1 are important, and the correlation function 
R(u),t) can be calculated perturbatively in A^/uj. In this 
regime the LDOS correlations can be found within the 
standard diagrammatic technique by expanding in the 
diffusive modest. 

In the most interesting case < 1, a naive applica- 
tion of the perturbation theory in fails in the "bosonic" 
sector. Smallness of f2 in the last term in (flU)) is com- 
pensated by large values of Xb, and thus the term QXb 
in Hb cannot be considered as a small perturbation. An 
accurate construction of the perturbation theory in the 
bosonic sector at £1 <C 1 is a subject of Sec. IIIII 

Though physically relevant Q is imaginary, see 
Eq. (|13p , we will perform calculations assuming that O is 
a real positive number and make the analytic continua- 
tion to imaginary Q at the final stage of the calculation. 
With real f2, both Hb and H F become Hermitian oper- 
ators with a discrete spectrum, that allows to write the 
spectral decomposition of the expectation value (flU]) : 



(^olXm^fe) 



(Xm|Xm)(0fc|0fe) 



^-2(E m +E k )t 



(21) 



where Xm(X F ) and ^fe(As) are the eigenfunctions of the 
Hamiltonians H F and Hb'- 



H F (j) m (X F ) = E m (f> m (X F ), 
H B <f> k (X B ) = E k (f> k (X B ). 



(22) 
(23) 
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In the next sections we perturbatively construct the 
eigenf unctions of the Hamiltonians p^|) and (j2"0)) and cal- 
culate the matrix elements 



(/o(g)|Xm){p^l(p)|^k) 

+ (qh(q)\Xm)(K (p)\<p k ) (24) 



(*o|Xm<M 



needed to evaluate (|2Ijl . 

In the fermionic sector, the perturbation theory is 
straightforward, while its application in the bosonic sec- 
tor is much more involved due to its non-compactness. 
Nevertheless we start with the bosonic sector in Sec. IIIII 
and use some notations introduced there in dealing with 
the fermionic sector in Sec. IIV1 



III. EIGENSYSTEM OF THE "BOSONIC" 
HAMILTONIAN H B AT SMALL Q < 1 

The main technical problem in studying the noncom- 
pact bosonic sector is that it cannot be decomposed into 
a "regular part" and a perturbation such that (i) the 
regular part is exactly solvable and (ii) the perturbation 
is small on the whole interval 1 < A < oo (in this sec- 
tion we denote the bosonic variable Xb by A for brevity). 
Therefore to have a well-controlled perturbation theory 
one should use different decompositions of the Hamilto- 
nian (fT9|) at small and large A in order to capture the 
behavior of the wave function near the singular points 
A = 1 and A = oo. 

At small A <C 1/fi, we will decompose the Hamiltonian 
as Hb = H B1 + Vbi with 



H B1 



I 



d x {X 2 -l)d Xl 



V B i = OA, 



(25) 



whereas at large A > 1 wc adopt the decomposition 
H B = H B2 + V B 2 with 



>(A- 



l?d> 



+ 0(A + 1), 

v B 2 = d x (\ + i)d x - n. 



(26) 



A. Small A: Legendre scattering states 

The Hamiltonian H B1 has a continuous spectrum 

k 2 1 



E 



(o) 



(27) 



(k is a real number) with the eigenfunctions finite at A 
1 given by 



4 0) (A)=P_ 1/2+lfc (A), 



(28) 



where P-i/2+ikW is the Legendre function (with such 
an index it is usually referred to as the conical func- 
tion). The standard perturbation theory with respect 



to the perturbation Vbi is ill-defined since the matrix el- 
ements of the perturbation diverge when integrated over 
the whole semiaxis A > 1. Nevertheless it is possible to 
construct the wave function perturbatively in the limit 
1 < A -c I/O. Using the properties of Legendre func- 
tions, one can show that Vbi acting on the bare function 
Lf, (A) shifts its "momentum" k by ±i: 



yB1 4°)(A ) = n (i * + ^ (A) + ( '*-^ I ft (A) 



2ik 



(29) 



Therefore we can search the exact solution as a formal 
series 



L k (X) 



Cn (n,fc)4°> (A). 



(30) 



Then acting by the Hamiltonian Hb and using that 
L^h n (X) is also an eigenfunction of H B1 we obtain a 
tridiagonal system of linear equations for the coefficients 
Cn(Q,k): 



1 



2 8 
fl(ik - n+ |) 
2(ik-n + l) 



— Ek 



Cn-l 



Q(ik-n - |) 
2(ik -n-1) 



c n+1 = 0.. (31) 



The solution (|30p can be multiplied by an arbitrary A- 
independent constant. To fix this freedom, we assume 



cb(n,A) = i 



(32) 



and will normalize the eigenfunction (f>k(X) later on, see 
Sec. UlTDl 

The system of equations (|3T|) allows to obtain the series 
expansion of c n (Q, k) in powers of 0, starting with fi' n '. 
In particular, 



ci(Q,k) 
c 2 (Q,k) 

C3 (n,k) 



in i{2k 2 + nk - 7)n 3 

2k ~ 32(k~i)k{k + i) 2 (k + 2i) 

{2k + 3i)n 2 
~ 16k{k + i) 2 + "" 
i(2k + 5i)n 3 



+ 



(33) 
(34) 
(35) 



96k(k + i) 2 (k + 2i) 

The coefficients with negative indices can be found from 

c- n (n,A)=c„(fi,-A0. (36) 

The energy Ek also acquires corrections in powers of Q. 
Using Eq. (|3Tj) , it can be expressed in terms of c±i(f2, k): 

co-, Q(ik + h) , , n(ift-i) , 



k 2{ik+ 1) v ' 2(ik - 1) 

Thus Ek is given by the series in il 2 : 

n 2 (5k 2 + 7)0 4 



E k - E k 



(0) 



4(/c 2 + I) 44(fc 2 + I) 3 (fc 2 + 4) 



(37) 
(38) 
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B. Large A: Bessel scattering states 

The eigenfunctions of the Hamiltonian H^ 2 which de- 
cay at A — > oo are expressed in terms of the modified 
Bessel function of the second kind (MacDonald function) : 



K 2l k{v) 



(39) 



where p = y/8Q(X + 1). The corresponding energy is 
given by Eq. (|27|) . In order to take the perturbation 

Vb2 into account we note that its action on b[ (A) is 
equivalent to shifting the index k by ±i: 



v B2 B { k °\\) = n 



(ik 



2ik 



(40) 

[Note the spectacular coincidence with Eq. ([29]) !] There- 
fore the solution of HsBfc(A) = EkBk(X) can be natu- 
rally written as a series 



B k (X) 



c n (Q,k)B^_.(X) 



(41) 



where the coefficients c n (f2, k) coincide exactly with the 
coefficients of expansion over Legendre functions in ([3H)l , 
and can be found from the set of linear equations (|31[) . 
The spectrum Ek as a series in f2 can be found from the 
same equation, see Eq. (|38l) . 



C. Wave functions of the discrete spectrum 

The Legendre and Bessel scattering states constructed 
above provide the solutions to the Hamiltonian (fT!?]) with 
the proper behavior at A — > 1 and A — > oo, respectively. 
These solutions exist for arbitrary "wave vectors" k. The 
wave function </>fe(A) satisfying both boundary conditions 
can be found only for a discrete set of k. Such a wave 
function can be conveniently represented as 



Lfe(A), A < A, 
CkBk(X), A > A* 



(42) 



with some arbitrary A» and a suitably chosen factor Cfe. 
Matching the left- and right-scattering states, Lfe(A) and 
CkBk(X), at A* provides the condition for the allowed 
values of k. 

In the limit fl <C 1, determination of the discrete spec- 
trum is simplified by the fact that there exists a broad 
intermediate region 1 <C A <C where both Lfc(A) 

and Bfc (A) can be represented by their asymptotic expres- 
sions. To find them we write each term in Eqs. (|30p and 
(I4ip as a superposition of right- and left-moving waves: 



4 0) (A) = 4 0)+ +4 l 



(0)- 



B?\X)=B^ + B 



(o)+ 



>(£>)- 



(43) 
(44) 



and represent the right-moving waves Lu and B^ + 
by the following hypergeometric series: 



L 



(o)+ 



~ (-i) m r(2^fc-m) (x + iy h - +lk ~ m 



!r 2 (± + ik- 



B 



(o)- 



m— 
°° ( _ 1 ^ m 

S^r r( ~ 2 ^" TO) [ 2 " (A + 1) 

m— 



(45) 
(46) 



The left-moving waves are expressed by the right-moving 

waves as = and B^ = B^ k + ■ 

The right-moving waves in the Legendre and Bessel 
scattering states should be matched by a suitably chosen 
factor: 



Lt(X)=C+B+(X). 



(47) 



We compute the value of C k in the following manner. 
First we substitute the asymptotic expansion Eq. (|45|) of 
the Legendre function into Eq. (j30|) . This yields a dou- 
ble series with the terms — (A + i)-i/ 2 +*-m-™. g^. 
stituting the Bessel asymptotics (14T>1) into the expan- 
sion (|4ip yields another double series with the terms 
— (A + i)- 1 /2+ifc+m-n T n or d cr to match these two 
asymptotic expansions, we pick the terms with n = —m 
from the Legendre series and with n = m from Bessel 
series. This way we obtain two different representations 
of the coefficient multiplying (A + i)^ 1 / 2 + l,c i n the com- 
mon asymptotic region 1 <C A <C 1/57. They are single 
series in c„, that is in powers of f2. By equating the two 
expressions we find C k . The other terms in the asymp- 
totic expansion of L k and B k will automatically match 
due to the properties of the coefficients c n (il,k). Thus 
we have 



C. 



4r(2ifc)(4ft) 1 / 2 " ifc 
T(-2ik)T 2 (l/2 + ik)' 

where the phase factor 



i*(fc) 



(48) 



c„(fi —k) 

_ n=0 



(-1)™ T(2ik- 



T(2ik) 



c n (tt, k) 



(-1)™ T{-2ik + n) 
~. T(-2ik) 



(49) 



has a regular series expansion in powers of fl 2 : 

i*(fe) = -, , %kn2 , 
4(fc 2 + l) 2 



(50) 



The waves in the opposite direction, L fc (A) and B k (A), 
are matched by the coefficient 



(51) 



The full function <fik can be matched only if the two co- 
efficients coincide, C k = C k , which is a condition for 
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the discrete spectrum. It can be conveniently formu- 
lated as the condition S(k) = 1 for the scattering matrix 
S(k) = C+/C^: 



where 



So (*0 - (40) 



-2ik 



Y(2ik)Y(l/2-ik) 



T{-2ik)T{l/2 + ik) 



(52) 



(53) 



Equation S(k) = 1 determines the values of the wave 
vector k corresponding to the states of the discrete spec- 
trum. For the low-lying states (k < 1), one finds 



irN k 



ln(l/4f2) : 



(54) 



with Nk = 1,2,... labelling the energy level (Nk = 1 
refers to the ground state). According to Eq. (fS"4|) . the 
lowest discrete fc's are nearly equidistant. This can be 
naturally seen if we rewrite the Hamiltonian Hg in terms 
of the angular variable 9 = arccosh A. The motion in the 
resulting potential will resemble the motion in a box of 
size ln(l/4f2), leading to the quantization rule (|54|l . 
At the points of the discrete spectrum, the two coeffi- 
coincide, and we denote them simply 



cients and C k 
by Cfc. A calculation of C\ = C^C^ 



using Eq. (|4*5|) gives 



Ck = Ok cosh7rfc, <7fc = (-1) 



(55) 



where Nk is the number of the discrete energy level, see 
Eq. 454}. 

Remarkably, the expression (|55|) does not have any cor- 
rections in fl and is exact. 



D. Normalization of (j)k(X) 

The eigenfunction ^fe(A) perturbatively constructed 
above does not have the unit norm (4>k\4>k)- Instead, by 
definition, it is normalized in such a way that its expan- 
sion pop in Legendre functions comes with co(f2, k) = 1. 
Therefore we further need to compute the norm of <^>fe(A). 
This can be done by writing the function (f>k(X) in the 
form (142|) and calculating the normalization integrals 
with [Lfc(A)] 2 and [E>k(X)] 2 . This program is performed in 
Appendix [A] leading to a surprisingly simple expression: 



(4>k\<Pk) = -i 



coth7rfc d\nS(k) 



2tt 



dk 



OO 

E 

n— — ck 



k + in 



(56) 



where S(k) is the exact scattering matrix (|52")l . 

(Strictly speaking, Eq. (|56"|) has not been proven. It 
has been verified to the order 0(f2 6 ), see Appendix IA1 
for discussion. Nevertheless, relying on the anticipated 
duality between the bosonic and fermionic sectors of a 
supersymmetric theory, we believe that Eq. (|56p is exact.) 



E. Overlap with the zero mode 

The last ingredient required for the evaluation of the 
general expression ([2"Tj) is the overlap integral (^o\4>k) 
with the zero mode \I/o- According to Eq. (12"4")) . this 
involves two matrix elements from the bosonic sector: 
(K (p)\<j)k) and (pK 1 (p)\(j)k} ■ They are calculated in Ap- 
pendix O and the result has the form: 



(K (p)\(j) k ) = cr fc 7r ^2 



c«(0, k) 
2^/Ti cosh7rfc ' 



(57) 



/ K ( \u \ c n (n,k)[l + 4(k + tn) 2 } 
{pKi(P)\<Pk) = crfcTr 2^ ; — ; ■ (5») 



n— — oo 



4V0 cosh irk 



The overall sign factor a k = ±1 defined in Eq. (JS5J) will 
drop from the final expression (|2ip . 



IV. EIGENSYSTEM OF THE "FERMIONIC" 
HAMILTONIAN H F AT SMALL Q < 1 

In the fermionic sector, the standard perturbation the- 
ory in small Q can be easily developed by treating the 
O-dependent term in the Hamiltonian (|20|) as a pertur- 
bation. However, instead of doing these routine calcula- 
tions, one can notice a formal duality between the bosonic 
and fermionic sectors. By substituting Xb >— > Xf, one 
finds Hb i—* —Hf- Therefore, the fermionic cigenfunc- 
tions Xm(Xp) can be readily obtained from the Legendre- 
series expansion (|3U| in the bosonic sector by imposing 
the condition of regularity at Xf = — 1. The latter condi- 
tion restricts k to — z(m+l/2) with m = 0, 1,. . ., render- 
ing Legendre functions P_i/2+ifc(A_p) be Legendre poly- 
nomials -Pto(Af)- 

Thus the fermionic eigenfunctions can be written with 
the help of the coefficients c n (Vl : k) as 

oo 

X m{X F )= c n {n,i(m+l/2))P m+n {X F ), (59) 



where we use the symmetry relation (J36J) . Since 
P-i-i(x) = Pi(x), we can rewrite Eq. ([59]) as a sum of 
Legendre polynomials with nonnegative indices: 



(60) 



where the coefficients 7;(£!,to) arc defined as 



ji(Q, m) = Q_ m (0, i(m+l/2))+c_i_i_ m (0, i(m+l/2)). 

(61) 

The eigenenergy E m is given by — Ek [see Eq. f|3T[) ] with 
the substitution ik i— ► m+ 1/2: 

„ m(m + 1) Q(m + 1) ,„ , , .. 

S m = — > „ -4 ^ci(n,i(m+ 1/2)) 



2m + 3 



2m- 1 



c-i(n,t(m + l/2)).. (62) 
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The wave function (|60|) is normalized according to 



1=0 



(63) 



at making the analytic continuation to imaginary Q, it is 
desirable to get rid of the sum over discrete spectrum in 
the bosonic sector existing only for real ft. This can be 
done with the help of the identity 



The overlap integrals of Xm{^F) with Io(q) and qh(q) 
can be calculated with the help of the expansion (|60[) and 
the following exact expressions: 



(I (q)\Pl{\ F )) 



21 + 1 I 2l+1 (4y/U) 



(64a) 



( qh{q )\ Pl (\ F )) = \/ —Q— . (64b) 



E- = - E / gn*) 



2vr 



dk 



(65) 



V. BEHAVIOR OF R(ui, t) AT SMALL w AND AT 
ARBITRARY t 



where we utilized the fact that the scattering matrix 
S(k) = 1 at the points of discrete spectrum and employed 
the Poisson resummation formula. 



A. General expression 

Now we are in a position to calculate A(t) with the help 
of the spectral decomposition (f2"Tj) . Since we are aimed 



As a result, the factor d\nS(k)/dk appearing in 
Eq. (f65|) cancels din S(k)/dk in the normalization (|56p . 
and the general expression (|2"Tj) can be written in the 
following form: 



A(u>,t)=Y2A^(u,t), 



A (0) (w,f) = 47r 2 Re V / dfcfc 

°° /-QO 

A( n> °)(o;,t) =47r 2 Re^ / dfcfc 

™_n J — OO 



tanh 7rfc 
-2— M mfe (fi)e 

COSh 7TK 



-2(B fc +B m )t 



tanli7rfc 

-S ,n (fc)M mfe (fi)e 

COSh 7TK 



-2(E fc +E m )i 



(66a) 
(66b) 
(66c) 



In deriving (|66c[) we used the symmetry of the scattering matrix, S~ n (k) = S n (—k), which allowed us to extend the 
integral to the whole real axis. 

In Eqs. (f6"6"]) . the scattering matrix S(k) is defined in Eq. (fB"2")) . the energies E k and E m are given by Eqs. (|3"T)) and 
(I62[) . and the function jVf m fc(f2, fc) coming from the matrix elements in Eq. (|21[) has a regular expansion in powers of 
fl which can be calculated with the help of Eqs. (|56|) . (|57j) , fj8 j) , ffl - <f63]) and ([64)1 : 



M mfe (ft) = 



' OO OO 

E E 

\i=0 s=— oo 



7;(^,m)c s (^,fc) 
21 + 1 



/ + j + (fc + is) 2 



l + --(k + lS ) 2 



(67) 



4fi 



^21 + 1 

kI=0 



7 ?(0,m) £ 



fc + is 



](n,k) 



Equations (|66[) provide a general expression for the func- 
tion A(u>,t) in the limit of small u>. 

The term A^(ui,t) given by Eq. (|66b[) is an even func- 
tion of fl responsible for the decay of the LDOS correla- 
tions at the localization length scale, see Eq. ([70)) below. 

The terms A^ (u>, t) with n > are given by the inte- 



grals (|66c|) with the strongly oscillating function 5q (fc) cx 
Sq (k) cx n~ 2mk . Therefore they can be naturally calcu- 
lated by deforming the integration contour to reach the 
saddle point 

*W = ^f^, (68) 
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In the process of contour deformation one would pick up 
the contribution of the poles of the integrand in (|66cp 
located at k = i/2, i, 3i/2, ... The interplay between 
the contributions of the saddle point and of the poles is 
responsible for the behavior of the LDOS correlator at 
the Mott scale Lm given by Eq. (J9j) . 

In Figure [IJ we plot the spatial dependence of the 
LDOS correlation function R(lo,x) = 1 + A(u>,i) calcu- 
lated numerically with the help of Eq. (|66|) . where it is 
crucial to take the real part and use the fact that £1 is 
purely imaginary according to Eq. (|13p . Different curves 
on the graph correspond to various values of u> in the 
deeply localized region, lo <gi A^. The smaller is the ra- 
tio w/Ag, the faster do the sums in Eq. ([66]) converge. 

The qualitative behavior of the LDOS correlation func- 
tion in the quasi-lD geometry coincides with the behav- 
ior in the strict ID geometr y 21 ' 22 : On increasing x, the 
function R(lo,x) first decays at the localization length 
and then reaches the uncorrelated value R = 1 at the 
Mott scale Lm- 



B. Leading order in to 

In this section we evaluate Eq. l[66|) analytically in the 
limit of vanishing frequency, w/A^ — * 0. Then the local- 
ization length and the Mott scale are well separated and 
one can obtain simple expressions for R(co,x) describing 
its decay at x ~ £ and further growth near x ~ Lm- 

It can be easily seen that the coefficients 7z(f2, m) de- 
crease with increasing of the fermionic "momentum" to: 
M mk (fl) = Oifl 2 " 1 - 1 ). Therefore, only the ground state 
(to = 0) of the fermionic sector contributes to A{uj,t) in 
the limit lo — > 0: 



.U„.(<2)= (i+1) 



+ i/4) 2 + o{n). 



(69) 



To the leading order in small to one can also replace 
S{k)^S (k), E k ->4 0) , E, 



E 



(0) 



since corrections 



to these quantities start with il 2 . 

Consider first the contribution of A^(uj,t). Since 
Refi = 0, the "dangerous" term 1/40 in ([69| does not 
contribute and we may safely put lo = 0: 



3 2 f°° 
A (o) (0,t) = 4tt 2 — / kdk 
ot J Q 



tanh irk 

2 ' 

cosh irk 



-(fe 2 + l/4)t 



(70) 



The terms A^ n ' (lo, t) with n > can be evaluated by 
deforming the integration contour to pass through the 
saddle point (|68p and picking up the contributions of 
poles to be crossed. Since the residue at the pole ip/2 
is proportional to Q np ~ , only the term with n = 1 is 
relevant to the leading order in lo: 



-I 



dk k ^ So(k) 

, cosh irk 



1 



X {k 2 + l/ 4 )2 e -(fc 2 + l/4)t + ( ) 



The value of the integral (|7T|) depends on the relative 
position between the saddle point fc* 1 ' = i{pM + iir)/2t 
and the pole i/2. They nearly merge at t ~ tM, where 



t M = 2Rcln— = 21n— ^ 

Ail LO 



» 1 



(72) 



corresponds to the Mott scale L M given by Eq. ©. Thus, 
within the present technique, the appearance of the Mott 
scale is related to competition between the pole and the 
saddle point. 

The contribution of the saddle point fc^ 1 ^ calculated 
with exponential accuracy is 



A^(LO,t) 



exp 



(t - t M f 



(73) 



The contribution of the pole i/2 to A^{L0^t) remains 
finite in the limit lo —t 0: 



A^(0,t) 



pole at i/2 



-l. 



(74) 



The pole contribution ([74]) should be taken into account 

only for t < tM, since at larger t, Irufe* < 1/2 and the 
pole i/2 should not be crossed in deforming the integra- 
tion contour to the saddle point. 

Various asymptotic regions for the correlator R{lo, x) 
are considered below. 



1. Decay of correlations at t < ^m/2 

At spatial separations x smaller than the Mott scale, 
t < tM, the LDOS correlation function (in the leading 
order in lo) is determined by the interplay of A^°\0,t) 
given by Eq. (f70l) . and the two contributions to A^ (lo, t) 
[Eq. ([73D and ([71)] ■ The pole contribution flTJ]) cancels 
1 in Eq. ([7]). Therefore, R(lo,x) is determined by the 
competition between Eqs. (|70|) and (|73|) . At t < tjv//2, 
the former dominates and the LDOS correlation function 
reduces to A (0) (0,t) given by Eq. (|70"|) : 



R(0,x) = 4tt 2 ^ 



k dk ■ 



tanh irk 

5 ' 

cosh irk 



,-(fc 2 + l/4)t 



(75) 



This function decays at the localization length scale and 
describes the limiting curve in Fig. (TTJ) as lo — > 0. 

The formula (|75p is well known in the theory of lo- 
calization in one-dimensional systems. It describes the 
decay of a single wave function, both in the 1D 15 ' 16 i 17 
and Q1D 1 ' 18 ' 24 geometries. In the two limiting cases of 
small and large t one finds: 



R(0,x) 




— — — + -rzr + . . . , t < 1, 



t » 1. 



(76) 



16t 3 / 
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2. Intermediate asymptotics at tu /2 < t < tu 

The LDOS correlator decays according to Eq. (|75|) only 
for t < t M /2. For larger values of t, the saddle point 
contribution (|73p should be taken into account on the 
background of the exponentially small A^°\0,i). They 
become comparable at t ~ ^a//2, which results in a min- 
imum on the R(lj,x) curve, of the order of e _ * A/ / 8 ~ 
(w/A^) 1 / 4 . At larger distances (*m/2 < t < tj^), the 
LDOS correlator increases as 

exp (77) 



TABLE I: The poles of the integrand (|66c|l contributing to 
the coefficient ai(t) in front of the ui 2 term in Eq. (|84|l . The 
index n labels the power of the S-matrix, and the index m 
labels the fermionic state in the matrix-element block M m k- 



pole 


n 


m 


exponent 


i/2 


1 





1 




2 





1 




3 





1 




1 


1 


e~ 2t 


i 


1 





e 3t/4 


3i/2 


1 





e 2t 



3. Behavior at the Mott scale, t sa tjvf 

In the vicinity of the Mott scale, at \t — t^j | ~ \ftu, 
the saddle point /c* 1 ^ approaches i/2, and it becomes im- 
possible to treat the pole at k = i/2 independently of the 
saddle point. To evaluate (|7T|) in this region we note that 
its first derivative over t does not have a pole at k = i/2 
and thus can be easily found with the steepest descent 
method: 



dAM(u,t) 



1 



: exp 



(t - t M f 



Integrating (|75|) in the vicinity of t M , we get 



R(ui, x) 



1 1 



rf 



t — tM 

2\/£m 



(78) 



(79) 



Equation (|T9[) describes the crossover from R = to 
the uncorrelated value R = 1 at the Mott scale. This 
crossover is characterized by a sufficiently narrow width 
y/tM <C f,w- The formula (|79[) exactly coincides with the 
result for the strictly ID geometry obtained in Ref. [2ll 



4- Beyond the Mott scale, t 3> tM 

Finally, we address the behavior at scales much larger 
than the Mott length. Here R « 1, and one is interested 
in the small correction R(ui, x) — 1. 

In the region in < ( < if f , reduces to 

A^'(u,x) which is entirely determined by the saddle 



point k^\ and we get 



see Eq. (IMl) . The main correction to R(uj,x) 
from the ground state with k ~ 7r/m(l/40): 



1 comes 



= I + 87r 3 Re 



ifctanh7rfc Mofc(0)e 



-(fc 2 + l/4)t 



cosh 2 7rfc d\nS(k)/dk 



(81) 

The scattering matrix for such a small k is S(k) rj 
5 (fc) » e 2ifcln ( 1 /4f2)^ Using the mat rix element from Eq. 
(|69[) and performing analytic continuation, we obtain 



^,) = 1 __ sm ^___ 



^ 3 



x exp — - H 



t , <M 47T 2 (^-^ 2 )i 



,1/ 



(82) 



In the region t M <^t ^tj^, Eq. ([82]) simplifies to 



•'Af 

2tt 6 

it(o;, x) = 1 g— sm 



M 

8ttH 



= -t/4+t M /2-47r 2 t/*l 



(83) 



Oscillations of the difference R(ui,x) — 1 with x can 
be easily obtained in the opposite limit of large frequen- 
cies, u ^S> Aj, where localization effects are weak and 
the standard perturbation theory can be applied. In this 
regime, R{lo, x) — 1 simultaneously decays and oscillates 
with the characteristic scale given by the diffusive length 
\J D/lu. At present, we do not know any explanation of 
the fact that these oscillations persist well into the lo- 
calized region or of the physical origin of the oscillation 
period scaling as t 3 M . 



R{lu,x) = 1 - 



32i 7 /2 



exp 



(t - t M f 
At 



(80) 



breaks down at longest distances, t 3> 



Equation 

t 2 M . Here the contribution of the saddle points fci" ' with 
n > 1 become equally important, indicating that the re- 
placement (|65p of the sum over discrete levels by an inte- 
gral is not adequate. Indeed, at t ^> t 2 M the discrete en- 
ergy levels in the bosonic sector are completely resolved, 



C. u correction to R{ui,t) 

Since to the leading order in small to the result f75|) 
coincides with its ID analog, it is instructive to study 
the subleading corrections which would reveal the differ- 
ence with the strictly ID situation. Such a difference has 
already been seen^i in the limit of small spatial sepa- 
rations, i « 1, and now we calculate the correction to 
Eq. ([75)) at arbitrary t. 
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At sufficiently small t, the function A{u>,t) can be for- 
mally expanded in powers of w/ A £ : 



where A (o) (0,i) is given by Eq. (1701 



(84) 



1=0 



where the coefficients a;(i) are polynomials of ln(Ag/w) 
of order two. The series ([54")) is asymptotic: at t ~ 
tut/ {21 + 1) the /'the term in the sum acquires a non- 
analytic dependence on u> [the upturn (|79[) at the Mott 
scale is an example of such a behavior for 1 = 0]. 
To the leading order we have 



a o {t) = -l + A^{0,t), 



(85) 



The subleading coefficient ai{t) can be obtained from 
the general expression (|66p by adding the contributions 
from A(°)(w,t) and A( n> ^{uj,t). The latter are calcu- 
lated by deforming the integration contour to the upper 
half plane and picking the residues at the relevant poles 
specified in Table Q] The resulting expression has the 
form 



ai{t) 



(2£ - 3<) 2 + 17C 



24 



24i 829 
18 + 432 



71" 
1 



— dkk 



IT 

24 
tanh irk 
cosh 2 nk 



it 971-4 at/4 , 15 
512 ' 72 

(fc 2 + l/4) 2 (40fc 4 + 56fc 



(2£ - i) 2 3£ - t 7 

— H 1 

6 16 

2(fc 2 + 1/4) 3 



24 
7) 



7T 

24 



18(fc 2 
I 



I) 2 



3(fc 2 + 1) 



t e -(^+l/4)t ; (86) 



where 



C = ln(A^) -2 1 = t -f-2 1 , (87) 

and 7 = 0.577 ... is the Euler's constant. 

Equation ([56")) is valid for t < £m/3 and provides the 
exact distance dependence of the lo 2 correction to the 
main frequency-independent part i?(0, x) calculated in 
the previous section, see Eq. ([75)) . At t ~ t AI /3, the sad- 
dle point fc* 1 ^ goes below the pole 3i/2, annihilating its 
contribution at t > tuj^- At longer distances, the lead- 
ing frequency-dependent correction is due to the saddle 
point in A^'{u),t) and is no longer proportional to ui 2 . 



D. Summary 

In Table [Hi we collect the results for various contribu- 
tions to R{u),t) obtained in the previous Sections. Be- 
sides that, the decay at the localization length scale and 
the behavior at the Mott scale are given by Eqs. ([75]) and 
(179)) . respectively. 



BEHAVIOR OF R(u,t) AT SMALL t AND AT 
ARBITRARY w 



VI. 



Here we present an alternative approach to calculat- 
ing A{ui, t) which does not require the knowledge of the 
eigensystem of the Hamiltonian H. This approach ap- 
plies to the case of small t and arbitrary ui and amounts 
to evaluating the expansion of A{ui, t) in powers of t. 



We expand the evolution operator e 2Ht in the general 
expression (jT5J) for A{u>,i) in series over t: 



A{u>,t) = J2A n {u)t n , 



n=0 



where 



A n {u) 



2 



Re(* |(2 J ff)"|* > 



(88) 



(89) 



It is convenient to represent the matrix elements in 
the above expression in terms of the variables p and q 
defined in Eq. ([IB]). The Hamiltonian H, Eqs. (H)) and 
(1201) . acquire the following form in this representation 



-d pP {K 2 -p 2 )d p + -d q q{q 2 

p q 



n 2 )d q + P 2 



8 

(90) 

where k 2 = 16fl = —Aiuj/A^. This operator is defined in 
the region < q < n, p > k. 

The ground state wave function <3>o [see Eq. fH|l ] con- 
tains modified Bessel functions of the arguments p and 
q with indices and 1 only. The result of the action of 
H on such a function can be reduced to a combination 
of the same Bessel functions multiplied by some polyno- 
mials in p and q. We arrange the resulting expression in 
the form 



(2Jf)*|*o) = E q l V k ct\p^,K)LM)K k {p). (91) 

i,fe=0 

The matrices & n > are even polynomials of q, p, and k. 
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TABLE II: Summary of the results for R(lo,x) at small frequencies and arbitrary distances, t — x/£. 



Distance 



to- independent term 



tj-dependent term (tM = 21n(Aj/tj) 3> 1) 



tM 



tM 



<t <t M 



t M < t < ti 



t 2 M <ZLt 



7T 7 / 2 
16t 3 / 2 

7T ' 

16i 3 / : 



exp 



• exp 



24A 2 



(t M - 3f) 2 exp(2i) 



~ exp 

_9/2,3 

" 32i 7 / 2 



■ exp 



4i 



2t+ 

a/ 



87T tj\,f t 



exp 



4i 

tM 



4 + 2 



4^(t 2 M -^ 
(*m+^ 2 ) 2 , 



Their first entries are: 



,(i) _ P 2 + q 2 




(92) 



To make sure that the approach developed in this Sec- 
tion is completely equivalent to the one adopted in the 
rest of the paper, we have calculated the first 15 coef- 
ficients A n (0) and verified that they coincide with the 
coefficients obtained by expanding Eq. (|70[) in a series in 
t. 



,(2) _ 



~2(p 2 - q 2 ) 
\{p 2 + q 2 )-2n 2 



(P 2 



2k 2 



To find A n (u)), we have to multiply Eq. (|9"Tj) by 'J'o and 
integrate with <1\f = qdq/AQ, and d\s = p dp/ Ail. The 
integrals of Bessel functions arising in this calculation 
reduce to the products of modified Bessel functions and 
polynomials, owing to the fact that the matrix c'™' con- 
tains only even powers of p and q. Thus we can calculate 
the functions A n (w) up to any n recursively. The first 
three terms of the expansion ([55)) are Aq(u>) = A(u>,0) 
given by Eq. (JSJl, A\(uj) — B(u>) given by Eq. ©, and 



A 2 (uj) = ^ + -^ Re [(3k 4 
3 15 L 



11k 2 + 16)(I 2 K 2 +I 2 K 2 ] 



An(n 2 - A)(I X K - I Kj)(I K + hKJ 
3(k 4 - 3k 2 + A)l\K\ - k 2 {3k 2 - 8)1$ K$ 



(93) 



where we omit the argument k of the Bessel functions for 
brevity. 

In the limit 10 <C Aj, we can expand the cumbersome 
exact expressions for A n (u>) in powers of u>. The first two 
terms of this expansion are [C is defined in Eq. (|87p ]: 



Ao(u) = - 

A 1 (uj) = - 
s 

A 2 (uj) 



36A 2 
2u; 2 
9A| 



(12£ 2 + 52£ + 43-3tt 2 ), (94) 
£ + 1), (95) 



15 



18Af 



(l2£ 2 + 40/3 + 41 - 3tt 2 ). (96) 



The same terms naturally appear in the expansion of 
R(u,t), Eqs. (|55"|) and (|55|) . in powers of small t. 



VII. DISCUSSION 

In this work, we have developed a general approach to 
calculating the correlation function of the LDOS in Q1D 
disordered wires in the unitary symmetry class using the 
quantum-mechanical reformulation of the supersymmet- 
ric cr-model. Our approach is convenient both for numer- 
ical evaluation of the correlation function (see Fig. [1} and 
for analytical treatment of various asymptotic regions. 

To the leading order in small u, we demonstrate that 
the decay of R(lo,x) at the localization length and its 
jump to the uncorrelated value of 1 at the Mott scale 
are described by exactly the same formulae in Q1D and 
strictly ID geometries (with the natural replacement 
£ - * £id and Lm — * -^m )■ However already the next- 
to-leading correction at small u> are differential w 2 ln 2 w 
in the Q1D geometry and w 2 lnw in the ID geometry. 

We believe that the difference between the subleading 
terms in the ID and Q1D problems should be attributed 
to the short-scale structure of the wave functions. In this 
case, the oj dependence of the subleading terms should 
be sensitive to the symmetry of the Q1D problem, with 
the orthogonal case being different from the unitary case 
considered in the present work. 

In the limit of very large spatial separations, x 3> 
Lm, we identify a number of supcr-Mott scales. At 
x ~ ^ln 2 (A^/w), the asymptotic behavior of R(uj,x) 
changes its functional dependence, leading to a faster de- 
crease of correlations. Even a larger super-Mott scale, 
^ln 3 (A^/a;), gives the period of decaying oscillations in 
R(u>, x) — 1, according to Eq. (J53J). The physical origin of 
these super-Mott scales still remains to be clarified. Note 
that the very-large- a; asymptotics of the LDOS correlator 
in the ID case has not been addressed in Ref. [U and it 
is thus impossible to make a comparison between the ID 
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and Q1D results in this region. 

The other interesting feature of the problem is the ap- 
pearance of the whole hierarchy of sub-Mott scales. Each 
time the saddle point crosses a pole, a narrow step 
of the width Sx ~ £ ln 1 ^ 2 (A^/o;) arises at some rational 
value of x/Lm- The erf behavior in Eq. (|79p. originating 
from the interplay between the saddle and the pole 
i/2 is the strongest of those steps appearing already in the 
leading order. The amplitudes of other, sub-Mott, steps 
are proportional to some powers of ui/A^ and therefore 
they can be seen only in subleading contributions. 

We believe that our technique developed for calculating 
the LDOS correlation functions can be used to find other 
low-frequency properties of Q1D systems. The long- 
standing problem in the field concerns the derivation of 
the low-frequency dissipative conductivity in Q1D wires. 
The main technical difficulty compared to the case of the 
LDOS correlations is that calculation of conductivity in- 
volves a much more complicated evolution operatori^ 
compared to (dHJ). 

Finally, it remains a challenging task to generalize our 
approach from the unitary to other symmetry classes. 
This is a sophisticated problem since the corresponding 
quantum-mechanical formulation involves more than two 
independent variables X 



the terms as 



where in 



ji.m— — oo 



/ - I B +I L 



(A2) 



(A3) 



we only keep zeroth-order [constants or logarithms] terms 
in (A* + 1) in the integrals 



C = C£Z» +1 J B^_ in (A)Bj ra (A) dX, ( A4) 

Inm = Za.+i J" 4+ m (A)4°] im (A) dX. (A5) 

There are two types of such terms: diagonal (n = m) and 
off-diagonal (n ^ m). They have different structure and 
are computed by different methods. 



1. Diagonal terms 

Diagonal terms give rise to contributions scaling as 
log SI. They are computed with the use of the normaliza- 
tion lemma in Appendix [Bl 

For the case of Bessel functions, the integral (IA4[) at 
n — m = is written as 
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APPENDIX A: NORMALIZATION OF THE 
EIGENFUNCTION <t> k (X) 

Here we compute the norm of the function 4>k{X) given 
by Eq. 02]): 



[L k {X)fdX + C 2 k / [B k (X)] 2 dX, (Al) 



where A* is an arbitrarily chosen value in the 
intermediate-asymptotics region 1 <C A* <C I/O. 

We can expand each of the two integrals in Eq. (|A1[) 
as a power series in (A* + 1), and their sum should be in- 
dependent of A* . So we can keep only zeroth-order terms 
in (A* + 1) from the start, to simplify the calculation. We 
denote the operation of singling out zeroth-order terms 
(constants and logarithms) by Z\*+i- 

To compute the normalization in a systematic way, we 
expand L k (X) and -Bfe(A) via (f3"0"| and (I4ip and combine 



Xj 

[B { ° ) ] 2 dX-Zy, +1 

1+e 



A* 



[B^fdX. (A6) 



We imply the limit e — > so that the regular part (with 
oscillating terms ommited) of the first term may be cal- 
culated using the normalization lemma in Appendix [Bj 



Reg/ [B^fdX 
_ T(2ik)T(-2ik) 

~ mi 



1 d Y(2ik) 

ln(2fie + — — In — f- 

v ' 2i dk T(-2ik) 



(A7) 

The second integral in (|A6[) is calculated using the ex- 
pansion of the Bessel functions (|4l)|) at f2(A + 1) <C 1. Wc 
again omit all oscillating terms and obtain 



Z A . +1 Reg/ [B h 

-l+e 



A* 



(0)l2 



dX 



_ T{2ik)T{-2ik) X* + l 

len ln ^- (A8) 

Adding the two contributions (|A7j) and (|A8j) cancels 
the dependence on e as it should be and gives 



r(2ife)r(-2ifc) 

mi 



(A9) 



13 



We repeat the same procedure for the Legendre part: 



2. Off-diagonal terms 



Z A .+i / [L^fdX 



A [4 0) ] 2 dA-Z A , +1 AiffdA, (A10) 
1 J\* 



with A — > oo. Applying the normalization lemma from 
Appendix |5J we find the regular (without oscillations) 
part of the first term: 



A 



Reg / [LffdX 



T(2ik)T(-2ik) 



1 1 

x 4 



T 2 (ik + l/2)T 2 (-ik + 1/2) 



, A + i i a ln r(2ifc)r 2 (-ifc + 1/2) 

n 2 + 2idfc U T(-2ik)T 2 (ik + 1/2) 



(All) 



In calculating the regularized second term of Eq. (|A10[) 
we use the expansion (T4"5j) and obtain 



Reg 



[4 0) ] 2 ^a 



r(2ijfc)r(-2ifc) 



r 2 (iA: + l/2)r 2 (-ifc + l/2) 
Adding the two terms, we find 

T(2ik)T(-2ik) 



■4 In- 



A + l 
A* + l' 



(A12) 



I L - 



r 2 (ifc + i/2)r 2 (-i/c + 1/2) 



x 4 



' A* + 1 1 d h T(2ik)T 2 {-ik+ 1/2) 

11 2 + 2idfc n r(-2i/fc)r 2 (ifc+ 1/2) 



(A13) 

Other diagonal elements and 7^ n may be obtained 
by analytically continuing Iq and Iq to fc — * k + in. 
Note that such an analytic continuation commutes with 
the Z A *+i operation. 

Finally, we add the Bessel and Legendre contributions 
together and, using the identity 



T(2ik)T(-2ik) 



T 2 {ik + l/2)T 2 (-ik + l/2) 
obtain 



coth(7rfc) 
Auk 



coth(Trfc) d 
Imik) = ^k-dk lnSoik) > 



(A14) 



(A15) 



where Sa(k) is the scattering matrix (|53[) without correc- 
tions in f2. 

The other diagonal terms are given by the analytic 
continuation 

T ,i\ t n \ coth(7rfc) 3 
Inn(k) = I 00 {k + in) = .. . — lnb (k + in). 



2ni(k + in) dk 



(A16) 



Off-diagonal terms form a series in with positive 
powers only. They involve overlaps of two functions 
which solve the same Hamiltonian at different (complex) 
energies. Such integrals may be conveniently computed 
by the Wronskian method, which results in 



(A* + l) 5 



(k + in) 2 — (k + im) 2 

dx-B^l n (A*)i?£ m (A*) - B^ in (X*)d x .B^ im (X*) 

(A17) 



and 



(A* + 1)(A* - 1) 
(k + in) 2 — (k + zm) 2 



-(") 



(A18) 



(0) 



dx*L<£l in (X 

These coefficients have the obvious symmetries 



(A19) 



By using the expansions (|45p. (|46|) of the Legendre 
and MacDonald functions, after some tedious combina- 
torics, one obtains the following form of the off-diagonal 
coefficients (we have used the relation (|A14[) to simplify 
products of Gamma functions): 



I B = (-1) 
which leads to 

'o. 



coth(7rfc) 



i\(k + i 



7l-\-m \ ' 



(A20) 



2 coth(?rfc) 



; n+ra \ ' 



m + n — 1 (mod 2), 
m + n = (mod 2). 



n\n — m\ (k + V- 2 , 

(A21) 

Thus the Legendre and Bessel contributions cancel 
each other at odd powers of SI and double each other 
at even powers of fl. 



Normalization 



The diagonal terms (|A16[) in the sum (|A2|1 contain log- 
arithmic derivatives of So(k+in) taken at different n. Wc 
may perform two transformation on those terms. First, 
with the help of gamma- function properties, we reduce 
the argument of So(k + in) to k acquiring new terms 
which form a series in f2. Second, we replace lnSo(fc) by 
\nS(k) according to Eq. (|52")l , again at the cost of adding 
a series in tt. Surprisingly, these two contributions ex- 
actly cancel the off-diagonal part of the sum (|A2[) ! Wc 
have observed this cancellation up to the order 0(il 6 ) 
and conjecture that it is exact at all orders. 
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The resulting expression (after canceling the off- 
diagonal terms) is proportional to din S(k)/dk with the 
full S-matrix l[S2jl : 



((f>k\<t>k) 



. coth7rfc d\nS(k) 



2tt 



Ok 



E ■ 

n— — oo 



in 



(A22) 



APPENDIX B: NORMALIZATION LEMMA 



and the other one by the basis of the discrete spectrum 
of the states constrained by the infinite wall at Xq: 

J^i^o^^Ux^kix') = 6{x-x'), x,x' > x Q . 

k 

(B7) 

First, note that in the case of an infinite wall, the nor- 
malization of the states of the discrete spectrum coincides 
with the regularized one: I Xo ,k = R.eglx ,k- Second, we 
may replace 



In one dimension, consider a particle moving freely in 
one asymptotic direction and fully reflecting of a poten- 
tial in the other direction. Then the normalization of 
the wave function in the region where the particle expe- 
riences reflection may be related to the scattering phase 
shift. 

Specifically, consider the one-dimensional Hamiltonian 



H 



P 

2m 



U{x) 



(Bl) 



with the condition on the potential U(x) — for x < 
and U(x) — > oo at x — ► oo. Then to any positive 
momentum k there corresponds a state of the continuous 
spectrum ^k(x) with the asymptotic behavior at x < 0: 



^ k {x < 0) = e 



ik(x-xo) _ —ik(x—xo)+i<Pz: (k) 



(B2) 



where xq is an arbitrarily chosen reference point, and 
<Px (k) is the scattering phase relative to this point. Ob- 
viously, (p xo (k) is a linear function of x with the slope 
2k. Consider now the normalization of the part of the 
wave function from xq to +oo: 



xo,k 



\V k (x)\ 2 dx. 



(B3) 



x 



Obviously, for xq < this normalization is a linear func- 
tion of xo plus a sinusoidal function. We define the reg- 
ularized normalization integral 



Reglro.fc = Reg 



\y k (x)\ 2 dx 



(B4) 



as I Xo ,k with the sinusoidal part subtracted. 
Then our normalization lemma states that 



R-eg/^fc 



dk ' 



(B5) 



To prove the lemma, consider the same problem with 
an infinite wall placed at a position xq < 0. Then, for 
any two points x,x' > xq we have two decompositions 
of the delta function: the one by the full basis of the 
continuous spectrum, 



dk 

—n(x)*k(x')=6(x-x'), (B6) 

Z7T 



E 



dk t^? 

27 E • 

71 = — CO 



^(k) dVx (k) 
dk 



(B8) 



and after subtracting (|B6|) from (|B7j) we obtain 
~dip Xo (k) 



dk 
2^ 



dk 



-(Reg/^)- 1 -1 



Ef dk 



for any x, x' > xq. (B9) 



Now we take the limit xq — > — oo. In this limit, the first 
term goes to zero as l/|a;o|, since both d(p Xo (k) / dk and 
KegI xok are linear functions of Xq. The second term goes 
to zero faster than l/|xo|, because of the rapidly oscillat- 
ing exponent. Therefore, the first term must identically 
vanish, which implies the result (|B5|) . 



APPENDIX C: EVALUATION OF {K \(pk) AND 
{pKM 

Let us start with the matrix element (Ko\4>k)- Using 
the perturbative expansions (|30|) and (j4Tj) of the wave 
function 4> k , we represent the matrix element in the form 



n=0 



M L kLWKo(p) 



+ C k d\Bl% n (\)K (p) 



(CI) 



with 1< A* < I/O. We then use the hypergeometric 
series (|45j) . (|46|) for Legendre and Bessel functions and 
integrate every term of these expansions with Kq (p) with 
the help of the following indefinite integral 



FF >yF ' 1 + a 1 \ ' 2 ' 2 '4 
^-iF 2 [l,— . (C2) 



(1 + a) 2 



The parameter k satisfies the quantization condition 
S(k) = 1, which ensures that the two expansions of the 
wave function in Eq. (|C1[) match at the point A*. This 
allows us to omit the dependence of the two integrals in 
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Eq. (|C1|) on A* by applying Za*+i operation denned in 
Appendix [A] 

Let us start with the Legendre part of Eq. (jCljl . We 



integrate every term of the series (j45|) for with i^o (p) 



using Eq. (|C2[) and expand the emerging hypergeometric 
functions. The contribution of the upper limit A* does 
not contain any constant or logrithmic terms in A* + 1 and 
therefore is totally annihilated by the Za*+i operation. 
At the lower limit A = 1, we have 



A* 



E 



dX L k 0) (X)K (p) 
T{2ik-m) (-l) m+1 (4^)' 



' T 2 a + ik-m + l) 1 

m,l—Q v ' 

4V^^i(4vn) 



f | — ifc + m — 
We first sum over the index m and obtain 

rA* 

dA L< 0) (A)Ko(p) 



2 J ftr (4\/«) 
— ik + m — I 

{k i ► -fc}. (C3) 



E 



i?r(2Z)!(4fi) 1 



^ sinh(27rfc)r 2 (i + ?:fc + /)L 2 (± -ik + l) 
j 2if (4\/ft) (2/ + l)4 v / 0-ffi(4\/f2)l 



[(i+«r+*»] j j 

+ {jfe h-> -fc}. (C4) 



Every term in the above sum over I is an odd function 
of k and cancels against its (k t— ► — fc) counterpart. This 
immediately yields zero: 



A* 



dA4 o) (A)if o (p) = 0. 



(C5) 



This identity also holds after the shift of the parameter 
t h k + in because the integral remains convergent in 
the limit A — > 1 and the operation Za*+i commutes with 
the imaginary momentum shift. 

Now we turn to the second integral in Eq. (|Cip and 
rewrite it in the form 



In the second integral we will omit the dependence on 
A* according to the Z>» +1 prescription. Substituting the 
expansion of the Bessel function (|46|) and using Eq. (|C2|) , 
we sec that every term of the hypergeometric series yields 
zero at the lower limit p = 0. Therefore the Bessel part 
of the matrix element is 



Z A , / dX B ( °\x)K (p) = 



16ficosh 2 (7rfc) 



(C8) 



This result holds for the shifted index k ^ k + in as well, 
despite the fact that both integrals in the right hand side 
of Eq. (|CJ6[) converge only when k is real. 

Collecting all the terms we obtain the matrix element 



mn cosh 2 (%k) ^ 



c„(fi, k). 



(C9) 



Calculation of (pKi\<f>k) is very similar. We separate 
the Legendre and Bessel parts of the matrix element as 
in Eq. (|Clj) . With the help of the result (|C5[) we again 
find that the Legendre integral yields zero 



A * 



W / d\Lf\\)pK 1 {p) 



= -2f>^Z A . + i^ dA4 0) (A)^ (p) = 0. (CIO) 



Next, we extend the Bessel integral to the whole axis, as 
we did in Eq. (|C6[) . and find 



Z A . +1 / dX B { °\X)pK 1 (p) = d\B { °>(X)pK 1 (p) 



(Cll) 



7T 2 (l+4fc 2 ) 

32ft cosh 2 (71-fc)' 



Za u j dXB { °\X)K (p) 



dX B^(X)K {p) 



Za* 



A* + l 



dXB^(X)K (p). (C6) 



The first integral (running from —1 to oo) is known to 
be 



dX B { °\\)K (p) 



160cosh J (7rfe) 



(C7) 



After the analytic continuation k i— > k + in and collecting 
all the terms, we get 



32iZcosh 7rfc ^— _ 

V / n— 



^c„(ft,£;)[l + 4(fc + z7i) 2 ]. 

(C12) 

Using the constant Ck from Eq. (|5"5")l we finally obtain 
the matrix elements in the form Eqs. (|57f and (|58|) . 
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